Matter-wave gap vortices in optical lattices 
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We predict the existence of spatially localized nontrivial topological states of a Bose-Einstein con- 
densate with repulsive atomic interactions confined by an optical lattice. These nonlinear localized 
states, matter-wave gap vortices, possess a vortex-like phase dislocation and exist in the gaps of the 
matter-wave bandgap spectrum due to the Bragg scattering. We discuss the structure, stability, 
and formation dynamics of the gap vortices in the case of two-dimensional optical lattices. 



Similarities between the physics of coherent light and 
matter waves can be successfully used to understand and 
predict nonlinear dynamics of weakly interacting Bose- 
Einstein condensates (BECs) in optical lattices [J, 0, 
0, In analogy with periodic photonic structures for 
light waves, such as photonic crystals, optical lattices 
form band-gap structures which modify diffraction and 
localization properties of BECs with both attractive and 
repulsive atomic interactions. In the case of a repul- 
sive BEC, optical lattices offer two intriguing possibili- 
ties. First, they enable nonlinear localization of a repul- 
sive condensate in the gaps of the Bloch-wave spectrum, 
without a confining harmonic potential 0, Q • Secondly, 
by manipulating the band gaps of the lattice Bloch-wave 
spectrum, the degree of BEC localization can be varied 
from a low density state spread out across the lattice to a 
condensate droplet tightly bound in a vicinity of a single 
lattice well. Both the BEC diffraction management in 
shallow one-dimensional (ID) optical lattices 0, 0] and 
the nonlinear localization of a BEC with repulsive atomic 
interactions |S| were observed in experiments. 

One of the fundamental questions is how distinctively 
superfluid properties of weakly interacting condensates, 
e.g., the dynamics of vortices p, are modified by an opti- 
cal lattice, and whether it is possible to spatially localize 
a BEC wavepacket with a nontrivial topological phase in 
a lattice. Several parallel studies in the physics of coher- 
ent optical jlflf and matter waves suggest that local- 
ization of topological defects is possible for an attractive 
condensate within a tightly binding lattice. The local- 
ized vortices in the attractive condensate would be anal- 
ogous to bright optical vortex solitons predicted to exist 
in total internal reflection gaps of dynamical photonic 
lattices [HI E| an d photonic crystal fibers with fo- 
cusing nonlinearity, and recently observed in experiments 
on photonic lattices [l4(. Localization of vortices in the 
repulsive BEC confined by a shallow lattice, i.e. when the 
standard mechanisms of spatial confinement are absent, 
remains an open problem. 

In this Letter we show that nonlinear localization of 
a repulsive BEC with a phase defect is possible within 
a complete Bragg-refiection gap of the Bloch-wave spec- 
trum. The resulting matter-wave gap vortices are spa- 



tially localized and dynamically stable. Due to their fi- 
nite spatial extent, which is much smaller than the di- 
mensions of the loose confining trap, the localized gap 
vortices can also be termed "bright" vortices, by analogy 
with vortex rings of light in nonlinear optical media. Sur- 
prisingly, the density and phase structure of a localized 
lattice vortex is determined by the features of delocalized 
"dark" vortices, which reside on a spatially extended pe- 
riodic background of matter-wave Bloch states. 

The framework of our analysis is set by the mean- 
field Gross-Pitaevskii model of an anisotropic BEC cloud 
loaded into a two-dimensional optical lattice potential, 
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This equation is obtained by assuming a tight confine- 
ment in the direction perpendicular to the lattice ("pan- 
cake" trapping geometry) and a standard dimensionality 
reduction procedure (see, e.g., Ref. J^)- It is made di- 
mensionless by using the characteristic length ql = d/ir, 
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of the lattice, where d is the lattice period, and m 
the mass of the trapped atoms. The wavefunction is 
scaled as ^ — > ^yJg^D, where the rescaled 2D inter- 
action strength is #2D = 4v // 2 / 7ra s /ai / , and a s is the s- 
wave scattering length. To eliminate any possible lo- 
calization effects which are not due to the optical lat- 
tice, we removed from consideration any additional trap- 
ping of the condensate in the directions of the lattice. 
The trapping potential is therefore taken in the form, 
V(x, y) = V [sm 2 (x + 9) + sin 2 (y + 0)}. 

Stationary states of a BEC are described by solutions 
of Eq. JU of the form: \P(r, £) = ip(r) exp(— ifit) : where 
\i is the chemical potential. Stationary states of the non- 
interacting BEC are found as -0(r) = 0k (r) exp(ikr), 
where the wavevector k belongs to the Brillouine zone 
(BZ) of the square lattice, and 0k (r) = 0k (r + d) is a 
periodic (Bloch) function with the periodicity of the lat- 
tice. The band-gap structure of spectrum, /i(k), of the 
atomic Bloch waves in the optical lattice is shown in the 
top panel of Fig. ^ m the reduced zone representation 
usually employed in the theory of crystalline solids. 

In a homogeneous BEC, a vortex with the topological 
charge one is characterized by a 2tt phase winding around 
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FIG. 1: Top: Matter- wave Bloch spectrum in a 2D optical 
lattice (Vb = 3.0) shown along the edges of the irreducible BZ, 
r(k x = o,k y = 0) -> = l,ky = 0) -> = 1,^ = 

1). Shaded and open areas show bands and gaps, respectively. 



Bottom: examples of the condensate wavefunctions, 



% (red 



color codes maximum density) at the marked points of the 
dispersion curves, corresponding to (a) off-site vortex in the 
middle of the first BZ, and (b) on-site vortex at the edge of 
the first BZ, shown together with their phase structures. 



the low density core on a constant density background. 
In a condensate loaded into an "infinite" optical lattice, 
a single-charge vortex can be stabilized in the form of a 
" dark" Bloch vortex on a periodic background of a Bloch 
state. Two examples of stationary Bloch vortices found 
by numerical solution 0,^1 of the time-independent ver- 
sion of Eq. (1) are shown in Fig. ^a,b). Bloch vortices 
on the background of the ground lattice state [i.e. at the 
lowest edge of the first band, Fig. [Ha)] have the winding 
structure of the phase typical for a charge one vortex in a 
homogeneous condensate, due to the trivial phase of the 
background Bloch state |l|,|5[. In contrast, the Bloch vor- 
tex residing on the excited Bloch state at the edge of the 
1-st BZ with nontrivial phase develops a highly unusual 
phase pattern [Fig. 1(b)]. The phase winding around the 
zero density centre is 2tt. 

The "dark" vortices on the ground (Bloch) state were 
found to be dynamically stable in [Tfl, [l7| , although their 
experimental observation could be difficult because their 
core size is of the order of the lattice period. Recent 
analysis of the Bose-Hubbard model for a strongly cor- 
related bosonic system with nearest neighbor repulsion 
in an optical lattice [HI also revealed the existence of a 
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FIG. 2: Variations of (a) normalized atom number and (b) 
peak density of the on-site (solid) and off-site (dashed) gap 
vortex within the complete gap (open area) at Vb = 3. Shaded 
areas are spectral bands. Examples of the spatial structure 
of the vortices at the marked points are shown in Fig. 3. 
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FIG. 3: (a,b) Density, \ip\ 2 , and phase structure of the off- 
site (a) and on-site (b) gap vortices at marked points on the 
existence curves in Fig. El (c) Typical profile of a strongly lo- 
calized off-site vortex (at \i — 4.0). (d) Example of a strongly 
localized broad off-site vortex (at ji = 4.0, n = 124), with the 
parameters off the scale of Fig. Eta), but peak density shown 
in Fig. Etb), point "d". Red color codes maximum density. 



single vortex with the particle density that can be either 
suppressed or enhanced in the vortex core. 

Quite apart from the spatially delocalized states within 
the Bloch bands described above, spatially localized states 
of BEC with repulsive interactions can only exist within 
the complete spectral gaps |2, l3- By imposing the 
vortex-like phase structure onto a spatially localized en- 
velope of a Bloch state at the M-edge of the spectral 
gap (i.e. in the vicinity of the point b in Fig. 1), we 
have found, numerically, different families of spatially lo- 
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calized "bright" matter- wave gap vortex solitons charac- 
terized by the dependence of the (normalized) number 
of atoms in the localized state, n = J \ip\ 2 dxdy, on the 
chemical potential [see Fig.EJa)]. Akin to the bright gap 
matter- wave solitons [HI , gap vortices have a clear atom 
number cut-off below which they undergo a delocalizing 
transition. The lowest energy families shown in Fig. [2| 
correspond to the on-site (localized on the lattice min- 
imum) and off-site (localized on the lattice maximum) 
stationary gap vortices shown in Figs.|Hfa,b) respectively. 
The peak density of the localized state is a function of 
chemical potential within the gap, as shown in Fig.[2fb), 
however the width (FWHM) of the high-density vortex 
core practically does not change across the gap. Gap vor- 
tices are strongly localized in the middle of the gap [see 
Fig. dfc)], but near the top (X) edge of the complete gap 
they develop extended tails that have the amplitude and 
phase structure of the background Bloch state. 

Strongly localized core of the gap vortex is similar to 
vortex " cells" that were predicted to exist in tightly bind- 
ing photonic lattices optically-induced in a focusing non- 
linear medium Here, we find that the "elemetary" 
vortex cells of two basic types form broader gap vortices 
by "tiling" the square lattice in a symmetric manner [c.f. 
Figs. 3 (c) and (d)]. The width of the gap vortex core is 
an additional parameter which characterizes higher-order 
families of the vortex states, however the peak density, 
n max = |^| max 5 °f the vortex core varies very little be- 
tween "elementary" and broad vortices [see point d in 
Fig. (2fb) indicating the peak density of the vortex in 
Fig. 3(d)]. The high-density core of a broad gap vor- 
tex preserves characteristic phase structure of a vortex 
"cell", with its 2ir winding structure, whereas its tails 
are characterized by a nontrivial phase pattern associ- 
ated with the nontrivial phase of the underlying Bloch 
state, as seen in Figs. 0^a,b). 

The localized states near the M edge of the complete 
gap can be described by the trial function: 

ij)(x,y) = Ar m exp(-r 2 /a 2 + i(/>) Bi(x,y), (2) 

where r 2 = x 2 + y 2 , and B\(x,y) is the Bloch state at 
the lower edge of the first gap. The Bloch state at the 
respective edge is always out-of-phase with the lattice 
potential, and due to separability of the lattice potential 
it can be well approximated as B\ « cos(a; + #) cos(y + 0) 
(see, e.g. 1]), where is a shift of the vortex center 
relative to the lattice minimum. The approximate form 
of the on-site and off-site vortex configurations is well 
captured by the trial function (|2j) with m = 2, 6 = 0, and 
m = 1, 6 = 7r/2, respectively. In both cases the necessary 
condition for the existence for a stationary state, i.e. zero 
linear momentum 20] J Im^^V^dxdy = 0, is satisfied 
in the limit of a ^> 1. The characteristic width of the gap 
vortex is given by a = pir, where integer p determines the 
size of the high-density vortex core in lattice periods. The 
examples of exact numerical solutions of well-localized 




FIG. 4: Snapshots of the density distribution for a dynami- 
cally unstable wide off-site vortex in Fig. 3(d) at evolution 
times t=252 (left) and t=600 (right). 



off-site vortices corresponding to p = 1 and p — 2 are 
presented in Figs. |3fc,d), respectively. 

By direct numerical simulations of Eq. (1) we have 
confirmed that both the off-site and on-site gap vor- 
tices can be dynamically stable. The phase distribution 
evolves and the phase dislocation precesses around the 
centre of the vortex, whereas the vortex magnitude and 
polarity of its charge are preserved in time. In contrast, 
broad gap vortices characterized by much higher atom 
numbers in the localized state are dynamically unstable, 
and experience rapid derealization at relatively long evo- 
lution times, as shown in Fig. 21 This process also de- 
stroys the phase structure of the gap vortex core. 

Finally, we discuss the feasibility of dynamical gener- 
ation of gap vortices in an experiment. To date, genera- 
tion of nonlinear localized states of BEC in the form of a 
bright gap soliton has been successfully achieved in a ID 
optical lattice near the edge of the first BZ. The key 
to the success of that experiment was a low atom num- 
ber (N ~ 10 2 ) that allowed generation of a fundamental 
gap soliton near the band edge. In contrast to the ID 
case, the existence of the particle number threshold for 
the localized states in 2D [ly means that the generation 
process can fail for insufficient initial atom numbers [2^ • 
Our analysis shows that the key ingredients to successful 
generation of the gap vortex are: (i) preparation of the 
BEC wavepacket at the edge M of the first BZ, (ii) ini- 
tial number of atoms well exceeding the threshold atom 
number for gap vortices, and (iii) the peak density of the 
BEC wavepacket being above the threshold peak density 
for a gap vortex. While the requirements (i) and (ii) are 
quite intuitive, the requirement (iii) is most stringent. 
To demonstrate the dramatic role of the peak density 
thresholds, we impose a 2tt phase winding ramp onto a 
broad wavepacket near the edge of the first BZ and ex- 
plore its evolution in the dynamical simulations of the 
Eq. (1). We therefore assume that, initially, the wave 
function has a form given by Eq. (J2J) with m, 6 = 0, 
and a typical density distribution shown in Fig. 5(a). 
The numerical simulations of the mean-field model in |2l[ 
suggest that the preparation (and localization) of the ini- 
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FIG. 5: Top: Evolution of the peak density of a BEC 
wavepacket with the initial peak density (circles) above 
threshold for generation of the gap vortices, (diamonds) be- 
low the gap vortex threshold but above the bright gap soliton 
threshold, and (squares) below bright gap soliton threshold 
(Vb = 3.0). Bottom: Snapshots of the condensate density at 
marked points: (a) an initial wavepacket, (b) generated on- 
site gap vortex with the phase structure (not shown) of Fig. 
3(b), (c) generated weakly bound steady state of four bright 
gap solitons, and (d) discrete diffraction. 



tial BEC wavepacket in the above form can be achieved 
via an adiabatic process of driving condensate to the M 
edge of the BZ with subsequent expansion in the lattice. 
The phase dislocation can be subsequently imprinted by 
a laser field |22|. The results of the numerical simula- 
tions for a shallow lattice (Vb = 3.0), with a = 4.57T and 
A = 1.2, 0.7, and 0.6 are presented in top panel of Fig. 
5 by curves (b), (c), and (d), respectively. In the cases 
(b), (c), and (d) the initial number of atoms exceeds the 
threshold atom number value for generation of an on-site 
gap vortex by 9, 2, and 1.7 times, respectively. However, 
the initial peak density in the case (b) lies above the gap 
vortex threshold, in the case (c) it is in the intermediate 
domain (shaded area in Fig. 5, top) between the thresh- 
olds for a gap vortex and a bright gap soliton, and in the 
case (d) the peak density is only marginally below the 
threshold for generation of a bright gap soliton. 

The results of the initial state evolution are shown in 
Fig. 5, bottom. In the case (b), the on-site gap vortex is 
generated in the transient process where the BEC den- 
sity distribution is evolving around a vortex state with 
the peak density of n max « 2.2, marked by a dashed 
line in Fig. 5. The one-to-one correspondence between 
the peak density of a stationary gap vortex state and its 
chemical potential [see Fig. 2(b)] allows us to determine 
the approximate chemical potential (with standard devi- 
ation) of the generated on-site vortex state [pointed at by 
an arrow in Fig. 2(b)], and identify it as a well-localized 
gap state. Only approximately 40% of the initial atom 
numbers are retained in the localized state. Growth in 
the initial number of atoms, i.e. due to a greater width of 
the initial wavepacket, leads to the increasing deviations 



from the average state, and to larger uncertainties in both 
/i and n max of the generated vortex. In the case shown 
in Fig. 5(c) a weakly bound state of bright gap vortices 
forms, with an arbitrary phase relationship between the 
soliton peaks. The retention of atoms in the localized 
states is almost 90%. In the case (d) the localized states 
do not form and instead the discrete diffraction of the 
initial state is observed. 

In conclusion, we have predicted novel types of local- 
ized vortex states of repulsive BECs in optical lattices, 
and identified the key requirements for their generation. 
Our results provide the next step in understanding the 
complex interplay between the supernuid behavior of a 
BEC with its inherent signatures, such as vorticity, and 
the nonlinear behavior of coherent matter waves in pe- 
riodic potentials. The experimental observation of pre- 
dicted structures seems feasible in view of the recent ob- 
servation of a 2D BEC in an optical surface trap, and 
the fascinating prospects of the creation of optical sur- 
face lattices |23| . Our predictions are also relevant to the 
light beams carrying phase dislocations in 2D photonic 
structures with a defocusing nonlinearity. 
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